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ABSTRACT 

In this Letter we present a theoretical scenario to explain the steep correlation between disk accretion 
rates and stellar masses observed in pre-main sequence stars. We show that the correlations and spread 
observed in the two best studied regions, p Oph and Taurus, can be reproduced by a simple model 
of single star formation from a rotating collapsing core and the viscous evolution of the circumstellar 
disk. In this model, the rate of rotation of the parent core sets the value of the 'centrifugal radius' 
within which the infalling matter is loaded onto the surface of the disk. As a consequence, the disk 
accretion rate measured long after the dispersal of the parental core bears the imprint of the initial 
conditions of star formation. The observed trend results naturally if, at the onset of the collapse, 
cores of all masses rotate with the same distribution of angular velocities, measured in units of the 
break-up rotation rate. 

Subject headings: accretion disks — (ISM:) dust — (stars:) formation 



1. INTRODUCTION 

Young stars are surrounded by disks for a large frac- 
tion of their pre-main sequence evolution. The accretion 
of material from such disks controls the disk properties 
and its evolution with time. Measurements of the mass 
accretion rate onto the star (M acc ) were obtained for T 
Tauri stars (TTS) in Taurus several years ago, showing a 
large sprea d of M aec for stars of similar properties (Gull- 
bring et al. QUI) . 

The discovery of disks around very low mass stars and 
brown dwarfs (BDs) and the measurements of the ac- 
cretion rate in these disks has provided an unexpected 
result, namely that, underlying the known large spread, 
there is a steep dependence of M acc on the mass of the 
ce ntral o bject, roughly _as M acc oc M 1 - 8 (Muzerolle et 
al. 120031: Natta et al. (2004). This correlation only be- 
came clear when a sufficiently large interval of stellar 
mass M* (roughly two orders of magnitude) was accessi- 
ble to observations. 

The physical origin of this correlation is unknown. 
There has been a number of suggestions that it may be 
due to mass-dependent stellar properties, such as its X- 
ray luminosity, which may affect the disk physical con- 
ditions an d the angular momentum transfer (e.g., Muze- 
rolle et al. 120031: Moh antv et al. 120051: Natta et al. l2006a|h 
iPadoan et al. I l)2005j) pr oposed a scenario of Bond i -Hoyle 
accretion to explain it. I Alexander fc Armitaee I l)2006f) 
showed that the observations can be accounted for by 
disk viscous evolution if the ratio of the disk to stellar 
mass at the end of the core collapse is fixed, but the 
viscous time scale is proportional to 1/Af»; they then 
interpret the spread of M acc values as an age effect. 

In this Letter, we propose a different, somewhat more 
basic scenario. We follow the collapse of a rotating molec- 
ular cloud core and the simultaneous viscous evolution of 
the resulting disk. Assuming that cores of vastly differ- 
ent mass have similar rotation rates measured in units 



of their breakup rate, a trend of M acc oc j\/f i- 8±0 - 2 j s nat- 
urally reproduced. The width of the correlation can, in 
this very simple model, be traced back to a spread in 
core rotation rates. In this way, the M acc measured in the 
pre-main-sequence phase bears an imprint of the original 
core properti e s, as wa s already discussed qualitatively by 
INatta et al. IpOfia. 

2. DATA 

Disk mass accretion rates have been measured in many 
objects in star forming regions. In two regions, Taurus 
and p Oph, they are known for a large number of objects, 
ranging in mass from few solar mass t o brown dwarfs. 
The p Oph sample (Natta et al. I2006a|) includes almost 
all the Class II objects (i.e., objects with circumstellar 
disks) in the region, selected from their mid-IR excess; 
the accretion rate is measured from the luminosity of the 
infrared hydrogen recombination lines. The Taurus sam- 
ple is more heterogeneus, including classical T Tauri stars 
for which the mass accretion rate has been determined 
from the optical and UV veiling, as well as a number of 
optically selected BDs where it has been derived by fit- 
ting the Ha profiles with magnetospheric accretion mod- 
els (see Muzerolle et al. 200H and r eferen ces therein) or 
from the Ca IR lines (Mohanty et al.|2005). The presence 
of a disk in the BDs is not always known. A comparison 
of the results in p Oph and Taurus can be found in Natta 
et al. (2006a). The accretion rates turn out to be very 
similar in the two regions; the relation between M acc and 
M* has practically the same slope (M acc cx M + 18±0 - 2 ) and 
the range of values for any M* is comparable. The only 
difference appears in the BD regime, where Taurus seems 
to lack high accretors. Given the different selection crite- 
ria for BDs in the two regions (optical colors vs. infrared 
excess), and the uncertainty of the M acc estimates in the 
BDs, it is presently difficult to evaluate the significance 
of this. 
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3. MODEL 

Our model is quit e simi lar to the model described 
by Hueso &; Guillot fl2005t h enceforth HG05; see also 
Nakamoto & Nakagawa 11995(1 . It has previously been 
used in the context of a theoretical study o f the crys- 
tallinity of dust (Dullemond, Apai & Walch 1200^1 . In 
this section, we describe briefly the model for the collaps- 
ing cloud and disk evolution, and re fer for more details 
to Dullemond, Apai & Walch 1)2006(1 and HG05. 

3.1. Collapsing cloud 

We start with a core with temperature T, turbulent 
velocity dispersion Av, mass M COTC and solid-body rota- 
tion rate O (radian/second). We assume that the core 
forms a single star surrounded by a disk. For simplicity 
we a ssume the core to be a singular isothermal sphere 
(Shu[l977). After onset of collapse (at which time we set 
our c lock t o t = 0) the infall rate is M infa n = 0.975 c 3 /G 
(Srm ll977j) . where c s is the effective sound speed, which 
combines the thermal and turbulent pressure. In the 
Shu model used here, the infall rate stays steady until 
the mass M core has accreted onto the star+disk system, 
and then it abruptly stops. 

The original core has a solid-body rotation rate f2, 
which is slower than the break-up rotation rate f^bi-cak = 
yj GM COIe / -Rcorej where i? C ore is the radius of the core. 
According to the singular isothermal sphere model, this 
radius is i? CO rc = GM core /(2c 2 ), so the breakup rotation 
rate can be written as inbreak = 2y/2c^ / 'GM COIC . 

Upon collapse the core matter will spin up due to the 
conservation of angular momentum, leading to the for- 
mation of a disk. The infalling matter falls onto the disk 
within the centrifugal radius R ce ntr(t), which increases 
with time as R ce ntr(t) oc t 3 until the end of the infall 
phase. In the following we will name the final (maxi- 
mum) -Rcentr(i) to be the centrifugal radius of the cloud. 
This radius is: 
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In our calculations, we will assume that the range of pos- 
sible values of u> is the same for cores of very different 
mass, or in other words, that the range in ratios of rota- 
tional energy over gravitational energy is independent of 
core mass. This means that small, low mass cores have 
lower values of f2 than big, high mass ones. 

An important consequence of this assumption is that, 
for a given value of u> larger mass stars form from cores 
with larger centrifugal radii. This will turn out to play 
an important role in the of the Af acc — M* relation that 
follows from the model, as we will discusse below. 

3.2. Disk formation and evolution 

As matter from the core accretes onto the disk, the 
disk itself evolves, expanding outwords. The surface den- 
sity £(-R, t) changes with time while matter continues 
to fall onto the disk within i? C entr(i) (see HG05). The 



disk evolution continues after the core has disappeared, 
until it has entirely accreted onto the star. The radial 
velocity vr is proportional to the viscosity coefficient 
v = akT m j fim p flK- Here a is the parameter of vis- 
cosity (Shakura & Sunyaev 1973), which we take to be 
constant, T m is the midplane temperature, and fix is 
the Keplerian frequency. The midplane temperature is 
determined by taking into account the heating due to 
viscosity and to irradiation by the central star and by 
the accretion shock. We include an effective viscosity 
caused by gravitational instabilities in regions where the 
Toomre parameter drops below unity. 

The evolution with time of the star+disk+core system 
is discussed in detail in HG05. At time to, when the 
core has accreted completely onto the star+disk system, 
the ratio of the disk to the stellar mass is an increas- 
ing function of i? cen tr; it exact behaviour depends on the 
interplay of the infalling matter with the viscous evolu- 
tion of the disk. If one forces i? C cntr to be independent 
of M core , then Mdi s k(io) oc M coro , as expected; if, as in 
eq.(2), i?ccntr oc M core , we find that the disk mass in- 
creases approximately as M*. The following disk evolu- 
tion at time S> t n is ap proximately self-similar (see, e.g., 
Hartmann et al. Il998tl . so that at any given time t^ to 
the accretion rate is to zero order proportional to A/disk- 

4. RESULTS 

4.1. Model parameters 

Since we wish to compare our model results to the Tau- 
rus and Ophiuchus samples we need to choose core prop- 
erties appropriate for these regions. The average mea- 
sured line widths l)Mvers fc BensoiJll983t iBelloche et alJ 
l2001|) lead to infall rates of M in m ^ 10- 5 M Q /yr (c s = 
0.35 km/s) in Taurus and Af infa n ~ 3 x 10 _5 M Q /yr 
(c s = 0.5 km/s) in Ophiuchus. We shall take these two 
values as representative for these two star formation re- 
gions. The ages of stars in Ophiuchus and Taurus show 
a peaked distribution with median values of ~ 0.5 Myr 
and ~ 1 Myr respectively, wit h a sm all tail of objects at 
older ages (Palla and Stahler 1200011 . With these typical 
ages and the mass infall rates chosen above, the star for- 
mation time scale is always much shorter than the age of 
the system. This is consistent with the lack of correlation 
between spectral type and stellar age in both regions. 

With c s chosen at some representative value, and M corc 
varied to cover a range of masses between 0.03 and 
2.5 Mq, the only remaining parameters of the model 
are a and lo. Rotation rates are measured in a limited 
number of cores, in general more massive than the cores 
considered here. The measurements are consistent with 
uo in the range 0.03-0.5, but there is a lar ge fra ction 
(~70%) of non- detections (Goodman et al. 1993; Bel- 
loche et al. 120011) . In this Letter, we will vary u> between 
a very low value (0.01) and a rather large one (0.3), and 
investigate the effect of this parameter on the M acc - M* 
correlation. 

Finally, we chose for the viscosity parameter a the 
value 0.01, which has been argued to fit well the proper- 
ties of TTS disks (Hartmann et al.[l998), but it remains 
a very uncertain parameter because the physics of the 
anomalous viscosity of disks is not yet fully understood. 

4.2. Results for the M &cc —M t relation 
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Fig. ^ compares the model results to the observations 
for Taurus (left panel) and p Oph (right panel). Each 
line corresponds to fixed uj and fixed Mj n f a n and time, 
as discussed before. For each u>, the models predict a 
relation between M acc and M* roughly as M^ 8 , which is 
in very good agreement with the trend in the data. Lower 
values of uj give lower values of M acc for fixed M», and 
we find that values between 0.3 and <0.01 can account 
for the spread of observed M acc at all Af* . 

The slope of the correlation Af acc -M„ is practically in- 
dependent of the exact values of any of the model pa- 
rameters. However, the value of M acc for any given M t 
depends on the value of the parameters; A spread in age, 
Afinfal] or a among objects in a region introduces a spread 
in the values of M acc for a given M*. From our model 
results, however, we estimate that the observed spread 
in M acc is likely dominated by the spread in uj. M acc is 
not a strong function of M; n f a n or a (less than linear) . It 
decreases significantly with time (rougly as i 1 - 7-1 - 8 ), but 
the observed spread of more than two orders of magni- 
tude (in both regions) requires that the objects are ho- 
mogeneously distributed over a factor ~ 20 in age; this 
is much more than observed, since only few objects in 
each region have ages so dif ferent from the median val- 
ues (Palla and Stahler l2000j ). 

5. DISCUSSION 

The models presented in this Letter account very well 
for the observed correlation between M acc and M* in 
both star- forming regions. Note that the model parame- 
ters have not been adjusted to fit the data, but are chosen 
a priori from independent observational evidence. 

The general agreement with the slope M acc oc M*- 8 is 
recovered for all values of the infall rate, age, etc., as 
long as uj is the same for cores of all masses. This is a 
direct consequence of the dependence of the centrifugal 
radius on M core (eq.5), as the accretion rate in the disk 
depends on the mass reservoir in the outer disk. If we 
assume i? ccn tr independent of M coro , we recover the well 
known result that M acc is roughly proportional to M* . 

The scenario we propose has a number of implications, 
which can be tested. Firstly, we assume that the same 
formation process holds for objects differing in mass by 
two orders of magnitude, from intermediate-mass stars 
to BDs. If so, cores capable of forming very low mass 
objects must exist in star forming regions, with infall rate 
and uj on average similar to that of more massive cores. 
Our conclusion that the spread of the measured M acc is 
due to a comparable spread of the angular momentum at 
the onset of the collapse and that cores of very different 
mass may rotate at the same fraction of the break-up 
speed can also be tested by future observations. 

Our models predict that the disk mass must depend 
strongly on the mass of the central object (roughly as 

M 2 ) and should be tightly correlated with M acc . A sig- 
nificant number of disk mass measurements exist only for 
the Taurus sample and are shown in Fig. [21 The y have 
been derived for TTS by Andrews and Williams (2005) 
by fitting sub-mm and mm fluxes with simple disk mod- 
els assuming a normalization of the dust opacity of 10 
cm 2 g _1 at 300 /im, and a gas-to-dust mass ratio of 100. 
For BDs, we show the recent measurements of Scholz et 



al. ((2006) , normalized to the same value of the opacity. 
While the distribution of the disk masses as a function of 
stellar mass does not conflict with the model predictions, 
the expected correlation with M acc is not seen in the 
data. Uncertainties in c s and age will likely spread our 
models, but will not undo the general predicted trend. 
Therefore there appears to be a disagreement between 
the models and the observations in this respect. This 
may be an argument against our model. However, mea- 
surements of disk masses are severely affecte d by ou r ig- 
norance of the dust properties (Natta et al. I2006bj) . of 
the gas-to-dust ratio and of their dependence on disk 
and stellar properties. In addition, in the BD case one 
should note that the current millimeter detections are 
still very uncertain and probably re presen t the outliers 
of the real distribution (Scholz et al. 12006(1 . 

The very simple models used in this Letter neglect sev- 
eral important aspects of the star formation process, in 
particular the formation from a collapsing core of multi- 
ple systems, rather than single stars. A survey of disk ac- 
cretion rates of binary TTS in Taurus (White and Ghez 
1200 ID shows that M acc of the primaries does not differ 
from that of single stars; Fig. ^includes a number of bi- 
naries, in which both components are accreting, and we 
find that for both primary and secondary M acc is well 
within the range of single stars with the same mass. In 
the context of our models, this is possible if the original 
core breaks into fragments which evolve independently, 
with a distribution in uj similar to that of the cores from 
which single stars form. However, we point out that the 
majority of the objects in Fig. ^ are single, and for them 
the core-collapse models we use should be appropriate. 

The similarity between the M acc -M» properties of Tau- 
rus and Ophiuchus originates, in our model, from a com- 
pensation between older age and lower Mi n f a n for Taurus. 
To disentangle these effects, one needs to measure M acc 
in older star forming regions, where the age effect may 
become important. This will be challenging, as fewer 
stars should be accreting at a detectable rate, and very 
low mass objects will very likely be below the detection 
limit. However, checking if the upper envelope of M acc 
is consistent with our predictions will be important. 

The model we propose explains the properties and evo- 
lution of disk accretion entirely as a result of the initial 
condition for star formation and the disk viscous evolu- 
tion. We neglect other processes, such as intermittent ac- 
cretion, photoevaporation, variability of the X-ray emis- 
sion from the central star, etc., which may affect at some 
level the snap-shot measurements of Af acc , but we claim 
that they do not control the basic observed trends. 

6. CONCLUSION 

In this Letter, we propose that the steep dependence 
of the mass accretion rate on the mass of the central ob- 
ject, seen in two star forming regions of different age and 
properties (Taurus and p Oph), results from the imprint 
of the initial conditions for the formation of the star+disk 
system. A simple model of disk formation and evolution 
from collapsing cores naturally yields the M acc oc M*- 8 
relation, provided that cores of all mass have a similar 
distribution of rotation rates in units of their break-up 
speed. We show also that varying uj within a range con- 
sistent with the (few) existing observations could account 
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Fig. 1. — Mass accretion rate M acc as function of the mass of the central object M* in Taurus (left panel) and p Oph(right panel). 
Symbols are measured quantities or upper limits. Each line is a model series for fixed u), as labelled, varying M cor c- In all model series 
Af; n f a n and the snap-shot time are kept constant (10 — 5 MQ/yr and 1 Myr in Taurus and 3 X 10 -5 M0/yr and 0.5 Myr in p Oph). 

reasonably well for the large spread of M acc observed 
at any M*, even though other factors may also cause 
a spread. The predicted correlation of M acc with disk 
mass is not well reproduced by the data; this requires 
further investigations into the disk model and the disk 
mass measurements. 
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Fig. 2. — Disk masses as function of M» (top panel) and of 
Aif acc (bottom panel) for the Taurus sample; binaries unresolved in 
the millimeter observations have not been included. Only a few of 
the Taurus BDs shown in the top Panel have measured M acc and 
appear in Fig. HI and in the bottom panel. Lines as in Fig. 1. Note 
that all models practically overlap in the bottom panel. 
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